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Introduction 

We  consider  the  straight  ahead  transport  equation  [l] 


EM 

5x 


+  ji(E)F  = 


-/i(E  )x 

K(t,E)F(x,t)dt  +  K(E0,E)e  0 


(1.1) 


F(0,E)  =  0, 

under  the  hypothesis  of  separability,  i.e. 
K(t,E)  =  G(t)H(E). 


We  assume  in  addition  that  the  interval  (0,E  )  can  be  partitioned 
into  a  finite  number  of  subintervals  (E^^  <  E  <  E^)  (n  =  0,1, ...,N) 
such  that  on  each  of  these  intervals  K(E,E)  may  be  closely  approxi¬ 
mated  by  some  constant  multiple  of  dpi/dE,  i.e. 


K(E,  E)  =  an  § 


(En  <  E  C  En-1). 


n- 


In  an  earlier  report  [jl"J  we  discussed  the  solution  of  Eq.  (1.1) 
under  the  following  circumstances; 

(i)  is  a  finite  integer  (n  =  1, ...,N+l)o 

(ii)  a  is  a  finite  integer  (n  =  1, . . . ,N+l,n^p),  a  =  oo, 

0  <  a1  <  • . .  <  ap+1  <  . . .  <  aN+1  <  0. 

(iii)  N  =  I,  a2  =  1,  arbitrary* 


3" 


(iv) 

(v) 


N  =  1,  =  -1,  arbitrary. 


ap  is  a  finite  integer  (n=ls...,N),  ~  ~  2.  * 


_  1 


0  ^  GC  ^  0  9  9  ® 


The  solutions  in  all  these  cases  are  expressible  in  terms  of 
polynomials,  exponentials,  Bessel  functions  and  confluent  hyper¬ 
geometric  functions  with  second  argument  2. 

We  shall  discuss  here  the  following  case, which  generalizes 
case  (v): 

(vi)  an  is  a  finite  integer  (n=l, . ..,N),  t 

0  <  oh  <  . . .  <  a  ,  ap+1  <  . . .  4  aN  4  0. 


When  p  ^  N,  the  solution  will  involve  the  function 

c . 

G(k,v)  =  e  u  IQ(ku)du 


(1.2) 


used  and  tabulated  by  Bice  [2J0 
2.  Properties  of  the  Function  G(k,v) 

Since  the  function  G(k,v)  is  relatively  new  in  analysis,  it 
is  advisable  to  record  some  of  its  properties.  Some  of  these  pro¬ 
perties  will  be  useful  in  the  next  section  where  we  actually  solve 
Eq.  (1„1)0 


Let  us  define 

H(k,v)  =  k  G(k,v), 

Then  it  is  true  that 
/V  r 


ume~u  I  (ku)du  = 


o 


k“+1 


<Om  H(k,v)-ve"v 


0 


The  proof  of  this  identity  can  be  made  quite  simply  by  a  mathe¬ 
matical  induction  on  m  if  one  replaces  u  by  t/k  and  multiplies  by 
before  applying  the  operator 

If  we  differentiate  Eq.  (1.2)  with  respect  to  k  and  integrate 
by  parts  twice  we  can  show  that  G(k,v)  satisfies  the  differential 


equation 


<k2-i)  || 


*kG  +  ve  j  k  Ic(kv)  +  I^(kv) 


(2.2) 


In  particular*  when  k  =  1  we  have  that 


G(l,v)  =  ve~v]l  (v)  +  I-,(v) 


(2.3) 


a  re 


suit  first  noticed  by  Bennett  [3]®  Since  it  is  obvious  that 


Gr(0,v)  =  1  -  e“v , 


it  follows  from  Eq.  (2.2)  that 


1-k  G  =  1  -  e 


-v  f  t  Io(tv)  + 


(2. It) 


(2.5) 


=  ve 


I0(tv)  +  I1(tv) 

V^t2 


If  we  let  v  approach  +00,  we  see  that 


G(k,  00)  =  l/Zl-k^ 


(2.6) 


when  k2  <  1. 

It  follows  from  repeated  differentiations  of  Eq.  (2.2)  that 
the  integral  on  the  left  of  Eq.  (2.1)  can  be  expressed  in  terms  of 
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G,  polynomials  in  v  and  k#  e~v  and  Bessel  functions  of  k¥  when 


k  <  1*  In  particular p 
¥ 


(1-k^)  ue~u  I  (ku)du  =  G(k,v)  -  ve~v|X0(kv)  +  kl1(k¥)r  (2*7) 


(l~k2)2  j  -u2e"u  Io(ku)du  =  (2+k2)G(k,v)  -  ve"v||2+k2+(l-k2)v|l0(kv) 

+  j^k  +  (l-k2)kvj  I-lOcv)  .  (2*8) 

The  computation  of  G(k,v)  for  small  values  of  k  can  be  per¬ 
formed  by  expanding  Io(ku)  into  a  power  series  and  integrating 
term  by  term.  In  this  manner  we  find  that  [2,  p.  !&) 


G(k,v) 


s.\2« 


2n  -u 
u  e  ' 


(nl  )£ 


=  (2>9) 


in  whicl 


A  (v)  =  1 
n 


1 +  v  +  £r + 


Since  0  <  A^(v)  ^  1  when  0  =  v  oo,  the  convergence  of  this  series 
is  dominated,  uniformly  in  v,  by  the  convergence  of  the  power  series 
for  (l-k2)“1//2» 

When  i-k2  «  1,  it  is  more  convenient  to  expand  I0(ku)  into 
a  power  series  in  1-k2  [>,  p.ite] 

CO  (-l)m(x-li?)mc^u)m]^(u) 

I_(ku)  =  y  — - — - — - -  > 

°  ml 

m=o 
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and  integrate  term  by  term*  If  we  observe  that 


e  Xm(u)umdn 


2m+l 


(this  can  be  proved  by  differentiating  both  sides  with  respect  to 
v),  we  then  see  that 


oo 

G(k,v)  -  ve*~Y  y 


+  V7) 


(2m+l)ml 


(2.10) 


If  we  make  use  of  the  inequalities  [te  p.  tel 


,(v)  ^  ev,  v*Vv  6 


m  -m 
m  e 


and  the  asymptotic  behavior  of  the  gamma  function,  we  can  show 
that  the  series  for  e"‘'vG(k,v)  converges  uniformly  in  v  and  k  when 
|k^-l|  =  1  -  6  <  1,  the  convergence  being  dominated  by  that  of  the 


series 


lk2-i!' 


Although  the  series  (2<>9)  and  (2.10)  both  converge  uniformly 
in  v,  neither  is  entirely  satisfactory  for  calculating  values  of  G 
when  v  is  large  and  k  is  near  one®  To  derive  an  approximate  for¬ 
mula  for  these  conditions  we  notice  that  as  a  consequence  of  the 
asymptotic  expansion  for  I  and  Eq®  (2*6)  we  have  that 


in  which 


(0,m)  = 


rcm+i) 

-  m) 


=  j/Ti-kJv 


^Oyn 

The  integral  with  integrand  e~v  t  can  be  reduced  after  m  inte¬ 
grations  by  parts  to  the  error  function  integral®  It  turns  out  that 
sufficient  accuracy  for  most  purposes  can  be  obtained  by  retaining 
only  two  terms  in  the  sum  and  so  we  shall  not  carry  out  this  reduc¬ 
tion  In  detail®  For  the  first  two  terms  we  have  that 


G(kj,v)~  — 


wkfl-kT 


.2 

e”u  dt  + 


te) 


-t2 
e  c  dt 


k  M-Vk: 


(5k-l)[l-H(t1)j 

i/01-0 


i/2rr v 


(2.11) 


in  which  H(t-)  is  the  error  function  extensively  tabulated  else¬ 


where  [  5  J 


H(t  )  =  ^ 

yur 


3 •  Solution  of  the  Transport  Equation . 

It  is  known  [ 1 J  that  if  f(s,E)  is  the  Laplace  transform  of 
F(x,E)j  then  under  the  hypothesis  of  separability. 


f(s.E) 


K(Eq,E) 


K(E« .E«  )dE* 
pTlPI+s 
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If  it  is  true  that 


KCE.E)  =  an  § 


(E  <  E  <  E  t). 
n  x  n-1  • 


then  it  is  true  that 


f(s,E)  = 


K(E0,E)(po+s) 

(/i+S)an+1 


“X-1  n-1 


vak+l~ak 


(E  <  E  <  E  , ) 
v  n  n-1 J 


Let  us  now  suppose  that  the  constants  an  are  integers  (n=l, . . .,N)S 
1 

that  =  -  ?r  ,  that  0  <  <  a2  <*..<'  and  that 

ap+l  ^  ap+2  ^  ° e  *  ^  ^  Oo  Then  if  E  >  E^,  f(s,E)  is  a  rational 

function  of  s,  and  so  its  inverse  Laplace  transform  can  be  readily 
computed  in  terms  of  exponentials  and  polynomials*  If  E  <  E,T,  then 


K(E  ,E)(n  +s) 


)tti 1  7T  (fk+s)ak+1  ak  (/% 

(/if  s )  ^  W+s )  ^ 2  (jy s )  VV* 


,ak+l“ak 


f (s.E)  = 


tyVS) 


K(E0,E)P(s) 


(ps)1'/2(pN+s)1/2(p  +s)aP  aP+1 


(3.1) 


in  which  P(s)  is  a  polynomial  in  s  of  degree  ap-cy^-l  since 
a.p-1,  a]£+3_“a]£  and  -Qpf  ar®  non-negative  integers.  The  inverse  La¬ 
place  transform  of  f  can  therefore  be  expressed  in  terms  of  the 
inverse  Laplace  transform  of 
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1 


(3.2) 


(p+s )  1/2"guN+s ) 1/2  (^p+s) 

and  its  derivatives  with  respect  to  pip.  The  inverse  Laplace 
transform  of  (3*2)  is  ^6,  No.  219,  No®  55!?J 


.x 


eV[//;V2(^)rIuHMar 


2.e 


-ppx 


^%"2Fp 


Fpjy  I  <FflrVcJ  • 


v* 

It  follows  that  the  solution  F(x,E)  can  be  expressed  in  terms  of 
the  function  G  and  its  derivatives  with  respect  to  k,  together 
with  exponentials*  Bessel  functions  and  polynomials.  In  view  of 
Eq.  (2.2),  it  is  not  necessary  to  use  derivatives  of  G  with  respect 
to  k. 


If  p  =  0,  i.e®  if  a1<  a2  <  ...  <  <  aN+;L  =  -  \  ,  the  above 

analysis  needs  to  be  modified  slightly.  In  this  case  we  have, 
when  E  <  Ejp  that 


f(s,E) 


K(E0,E)P(s) 

(pQ+s )  lc"aI(/u+s  ) 1/2  (jUjj+s  ) 1//2 


in  which  P(s)  is  a  polynomial  in  s  of  degree  -a^.  The  formal  mani¬ 
pulations  to  invert  f(s,E)  are  thus  identical  with  those  to  invert 
(3*1)  if  we  set  p  =  0,  a  =  1. 

We  shall  now  apply  this  general  analysis  to  a  particular  case 
in  which  the  details  are  relatively  simple. 
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Let  us  suppose  that  N  =  1,  p  =  0®  Then 


-a- 


f(s.E)  _  _ _ _ 

K(E  pE)  1»gu  -j  /p  *■? /p 

0  (p0+s)  (p+s)1^2(p1+s)1/ 


(p  o+s )  (ja+s ) 1^2  (p1+s)V2 


y-1  /-ttin  /Pi-Pq' 
Z__  rK-a^r}!  U  +s 
r=o  \ 


»a- 


r 


(-aiH(p0-pi)r  ar 


1/2. „  *=.-,1/2 


™  (rl)  (-a.j-r)l  <)uQ  I  (pQ+s)(p+s)  /  (p^+s) 


Therefore 


F(x.E) 
,S7 


(-a.^)  I  (p-Q~P-i)  3r 

F=“o  ^ r  1  )2  ^ ,”ai“r  ^ 1  K 


/u0(y^)  “  ^rn)jrT  fi 


\ 


^  2(P-Pl)y 


0 


dy 


nx 

e  0  F(x„E) 


1 


x„E )  _  V 

^  fe(rl)2(-arr)l 


x 


(y-x)~  e 


^P+Pl-2>*0)y  ' 


|(p-|i1)y  (  dy. 


0 


There  is  little  to  be  gained  by  using  Eqs.  (2.1)  and  (2.2) 
to  reduce  this  expression  in  general.  The  method  can  be  illustrated 
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in  case  a.^  = 


■  lo  Then 


11X 

P  „  FjXj,.! 

k(e",e T 


l-(p1-/i0)(y-x)  dy 


2l1+^r^o)x 

m-2F0 


5  *  5(PTr2Po)x 


(PrfV  [  ye  a(F  2|a°)yi0 '^-n^yjay. 


If  we  now  use  Eq0  (2<>7)  we  find  after  a  little  manipulation  that 


e  °  F(x.E)  _  J2_ _ 

K(e“E)  /i+p2“  2pc 


(1  -  +  2hv)G(k,v)  + 


2hve”v  [l0(kv)+kl1(kv)] 


in  which 


FrM0 

/i+pr2pio  * 


M-Mi 

V  —  .  '  TT 

Wr2Mo  ' 


£  (p.+fi1-2fiQ)xs 


k7  -  ^ 


This  result  has  already  been  applied  to  calculate  the  transmission 
of  2  Mev  Y-rays  in  tungsten  [_?]• 
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